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SUMMARY 

A dye-tracer technique has heen developed whereby the quantity of 
dyed water collected on a blotter-wrapped body exposed to an air stream 
containing a dyed-water spray cloud can be colorimetrically determined 
in order to obtain local collection efficiencies, total collection effi- 
ciency, and rearward extent of impingement on the body. In addition, a 
method has been developed whereby the impingement characteristics ob- 
tained experimentally for a body can be related to theoretical impinge- 
ment data for the same body in order to determine the droplet size dis- 
tribution of the impinging cloud. Several cylinders, a ribbon, and an 
aspirating device to measure cloud liquid-water content were used in the 
studies presented herein for the purpose of evaluating the dye -tracer 
technique. Although the experimental techniques used in the dye-tracer 
technique require careful control, the methods presented herein should be 
applicable for any wind tunnel provided the humidity of the air stream 
can be maintained near saturation. 


INTRODUCTION 

In the design and evaluation of thermal icing protection systems 
for aircraft, one of the more important factors in the determination of 
heating requirements and the chordwise extent of heated area is the 
cloud-droplet impingement characteristics of the various components 
requiring protection. Previous investigations (refs. 1 to 11). have 
generally been confined to analytical solutions of the droplet impinge- 
ment characteristics. These solutions involve the use of differential 
analyzers or mechanical analogs for computation of the droplet trajec- 
tories once the flow field about the component has been calculated. 
Bergrun determined a limited number of trajectories for an airfoil in 
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reference 4 using step-by-step calculations. In addition, in reference 
12 Bergrun derived an empirical method permitting estimation of the 
droplet impingement characteristics on arbitrary airfoil sections. This 
method, however, appears to be more applicable to Joukowski-type air- 
foils than to low-drag or high-speed airfoils because the basic data 
used in developing the method were obtained for four Joukowski airfoil 
sections, but for only one low-drag section. 

Because of the labor and time required to calculate the complex flow 
fields about some bodies and the droplet trajectories through these flow 
fields, it is often impractical to obtain droplet -impingement character- 
istics theoretically. It is therefore desirable to obtain an experi- 
mental method for the determination of droplet impingement characteris- 
tics and thereby avoid the calculation of complex flow fields. 

A series of investigations were initiated at the NACA Lewis labo- 
ratory and under NACA contract at the Daniel Guggenheim Airship Instir 
tute (Univ. Akron) to determine a suitable means for experimentally 
studying droplet impingement characteristics on a body at above-freezing 
air temperatures. Initial studies at the University of Akron were con- 
cerned with developing a technique to obtain a record or trace of droplet 
impingement on an arbitrary body. Subsequent studies at the Lewis lab- 
oratory were concerned with improvements in the technique and the corre- 
lation of experimental impingement results with theoretical studies. 
Because available instruments were inadequate to determine accurately 
droplet size distributions, a method was developed by which the droplet 
size and distribution in the cloud were ascertained from impingement 
records on several bodies for which the theoretical trajectories are 
available. 

This report presents a description of tracer techniques developed 
by the NACA to obtain the experimental impingement characteristics of 
arbitrary bodies. In addition, a method is presented for applying theo- 
retical impingement data (refs. 1 and 2) to the experimental impingement 
data from cylinders to determine the droplet size distribution of the 
impinging cloud. . A review of basic impingement relations precedes the 
main body of the report to clarify and provide a basis for the develop- 
ment of the methods herein described. 


BASIC IMPINGEMENT RELATIONS 

A body moving through a cloud of water droplets will not, in gen- 
eral, intercept all the droplets originally contained in the voliime of 
air swept out by the projected frontal area of the body. If the body 
is considered to be fixed in a reference system, then, at a large dis- 
tance ahead of the body, the droplets and the air both approach in 
straight line paths with a free-stream velocity Vq. ' (All symbols used 
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herein are defined in appendix A.) As the droplets and the air reach 
the vicinity of the hody, the air streamlines deviate from their orig- 
inal straight line paths and flow around the body. The, air flow field 
is dependent primarily on the shape of the body and the orientation of 
the body with respect to the direction of the free-stream velocity. A 
droplet, because of its momentum, tends to maintain its original straight 
line path toward the body; however, the drag force imposed by the rela- 
tive velocity of the air with respect to the droplet as the air stream- 
lines deviate tends to cause the droplets to follow the air streamline. 
The combined effect of the momentum and the drag determines the trajec- 
tory of the droplet. For the case of small droplet size, high air vis- 
cosity, low free-stream velocity, and high air density, the effect of 
the droplet drag will exceed that of the droplet momentum; consequently, 
the droplets will tend to follow the curved streamlines around the body. 
In this case, many of the droplets which were originally directly ahead 
of the frontal area of the body will miss the body entirely. For the 
condition of large droplets, low air viscosity, high free-stream velocity 
and low air density, the effect of the momentum is greater than that of 
the drag, and the droplets will tend to maintain their original straight- 
line trajectories. For this case, only a few of the droplets which were 
originally directly ahead of the, frontal area of the body will diverge 
enough from their straight-line paths to miss the body. 

The basic differential equations of motion for a droplet in an air 
stream approaching a body are presented in dimensionless form in refer- 
ences 1 to 3. For a body configuration where shape and orientation with 
respect to free-stream velocity direction are given but size is variable, 
the theoretical trajectory of a droplet is determined by any two of the 
four impingement parameters K, (p, Rsq, and \|f. These dimensionless 

parameters are interrelated (K<p = Roq^, = <p/K) and include such 

factors as droplet size, body size, air viscosity, droplet density, air 
density, and free-stream velocity (appendix A). 

Body of arbitrary shape in cloud of uniform drop size . - The basic 
definitions of the dependent impingement parameters - limit of impinge- 
ment, total collection efficiency, and local collection efficiency - can 
be expressed in general equations for bodies of arbitrary shape and 
orientation. A general case for impingement considerations is that of 
an arbitrary body placed in an incompressible air stream which contains 
droplets of uniform size. For simple bodies, the flow field may be ex- 
pressed mathematically, and the equations of droplet motion may be solved 
by the use of a mechanical analog or a differential analyzer. The direct 
resialt of these confutations is a trace of the droplet trajectory for 
each initial ordinate measured from a reference line at an infinite dis- 
tance ahead of the body. For llltistrative purposes, the droplet tra- 
jectories for a sinqjle two-dimensional body are shown in figure 1. The 
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dashed traces represent typical trajectories for a fixed set of impinge- 
ment conditions (described by K and <p). The solid lines are typical 
traces for another set of conditions (other values of K and <P ) for 
which the momentimi effect is greater and the droplet follows a straighter 
trajectory. For a given set of K and «p conditions, the droplets 
which start out at progressively larger initial ordinates (yQ's) will 

strike the body progressively farther back along the siirface until a 
maximum value y^ ^ is obtained. The trajectory starting at Jq ^ 

will be tangent to the surface of the body at point p, which is a dis- 
tance Sm (measiured along the surface) from a reference point on the 
body. All droplets which have initial ordinates greater than yo^m 

an infinite distance ahead of the body (x = - •) will miss the body 
completely. Therefore, Sj^ denotes the distance along the body s-urface 

to the limit of impingement or the point on the surface beyond which no 
droplets will impinge. For arbitreiry three-dimensional bodies the rear- 
ward limit of impingement will vary everywhere around the surface 
periphery. 

For an arbitrary three-dimensional body in a cloud of uniform drop- 
let size, the total collection efficiency Ej^ is defined as: 


(amount of water caught in a imit of time by a body) 
(amount of water which could be caught in that unit 
of time if the trajectories were all straight lines 
parallel to the free-stream velocity direction) 


( 1 ) 


or 



( 2 ) 


where A^ is the projected frontal area of the body and Aq^jj^ is the 

area through which all the droplets hitting the body must pass at some 
point a large distance ahead of the body, both areas being normal to 
the free-stream velocity direction. If the numerator of equation (l) 
is denoted as the total water collection rate Wj^ (ib/hr), the equation 

may be written as: 


= 0.329 Vq w^ Af (3) 

The local collection efficiency p for specified impingement con- 
ditions (K, <p) is defined for any point on the body sxirface as: 
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(amount of water caught on infinitesimal area of surface 

in a given unit of time) 

(amount of water which could he caught in that time on 
that area if the trajectories were straight lines paral- 
lel to the free-stream velocity direction and if the 
area were oriented to be perpendicular to the 
trajectories) 


or 


3 = 


lim 

AA ♦O 
s 


AAo 

AAg 


dAQ 


(5) 


where AAg is an increment of area on the surface of the body and AAq 
is the increment of Aq (which is the area perpendicular to the air 
stream) through which all the droplets impinging on AAg must pass. 

The value of 3 at any point on the body determines the local impinge- 
ment rate Wp, defined as the amount of water collected on a small area 

at that point per unit time. In terms of the local collection efficien- 
cy, ¥p is defined as: 


Wp = 0.329 Vq 3 


( 6 ) 


The 3 values over the surface are also, useful in some instances for 
the determination of the total collection efficiency Ej^. An integra- 
tion of the local collection efficiencies over the surface area subject 
to impingement and division by A^ yields the following relations: 


Af 


/ 


3 dAg = 


Af 


/ 


dAo 

dA, 


dA 



f 



(7) 


If equation (?) is multiplied by 0.329 Vq w^ kf, it is apparent from 
equations (3) and (6) that 

/ Wp dAg = (8) 

Effect of nonuniform droplet size distribution on impingement pa- 
rameters. - In many cases droplets of various sizes are intermingled in 
a cloud whether it is naturally or artificially produced. Various graph- 
ical representations may be employed to describe the distribution of 
droplet sizes in a cloudj these representations employ either (l) a block- 
type distribution or (2) a smoothly faired distribution c\irve. These 
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droplet size distributions are often shown as the variation of emulative 
liquid-water content w as a ttinction of droplet radius a. The dimen- 
sionless form of their plot is n = w/w^ as a function of- a/aj^g^. By- 

definition, the volume-median droplet size is that droplet size 

for which half the total liquid-water w.^. is contained in droplets 
larger than a-med half the total water is contained in droplets 

smaller than a^g^. 

A series of dimensionless hypothetical block distributions is pre- 
sented in reference 2; an example of such a distribution (Langmuir D- 
type) is shown in figure 2 together with a smooth curve approximation. 
Both block and faired curve droplet size distributions may be employed 
in the analysis of the Impingement of cloud droplets on various bodies; 
the choice is in most cases a matter of convenience; 

The dimensionless dependent impingement paramaters s^L, Ejjj, and 

P for a body in a cloud of nonuniform droplet , size must be treated in 
a manner somewhat different from those for a body in a cloud of uniform 
droplet size. The limit of impingement-, which can be denoted by Sj^/L 

for a single droplet size, must now be defined as the limit of impinge- 
ment caused by the maximum droplet size present in the cloud, and is 
therefore denoted as (Sjn/L)^gjj. The parameters denoted by P, Wp, E^, 

and Wj^ for the case of a single droplet size are weighted to account 

for the effect of all droplet sizes present in the, distribution. The 
weighted paramaters used with a droplet size distribution are denoted 
herein as p, Wp, E^, and Wj^. 

The process of weighting these parameters when the droplet distri- 
bution is specified and the theoretical dependent impingement parameters 
Sm/L, P, and Em are known as a function of the independent impingement 
parameters K and (p is an averaging process. This process defines 
Em and ^ so that the equations 

Wm = 0.329 Afl^VoW^ (9) 

and ' ' . 

Wp. =.,0.329 Vq p (10) 

are satisfied. Equations (9) and (lO) for a nonuniform droplet size 
distribution correspond to equations (3) and (6), respectively, for a 
distribution of uniform droplet size. 



NACA TN 3338 


7 


If a cloud of block-type distribution such as the solid lines in 
figure 2 (with aj^g^ and w^ specified) impinges upon a given body 

for which theoretical trajectory calculations are available, the weighted 
total collection efficiency Em may be defined as an average En,: 

^ (^,1 ^^1 + ^,2 ^^2 + • • • Em ,7 

X 

= ®m,l '^l + ^,2 ^^^2 + • • • Em,7 ^'^7 

Every droplet size group in figure 2 has a corresponding finite value of 
partial water content Aw (or Zin) and droplet size a. For given oper- 
ating conditions, every droplet size specifies a K value which, along 
with <p, determines E^ for that droplet group. These E^ values may 

be then used in equation (ll) to determine 11^. Equation (ll) may be 
graphically represented by a summation of the areas in figure 3(a). 

If a faired-type distribution such as the dashed line in figure 2 
is- specified, equation (ll) takes the form: 

r v=wt rn=i 

Em = Em dw = Em dn (12) 

^t Wv=0 Jn=0 

For this type distribution Em may be plotted as a function of w (or 
n) as in figure 3(b) and integrated graphically to obtain Em. Then 
Wm may be determined from equation (9). 

The weighting of local impingement rate parameters such as p or 
Wp is accomplished in the same manner as the weighting of Em except 

that the averaging process must take place at every point on the surface 
of the body where P is desired since p varies along the surface of 
the body. The quantity P is the average P at the particular loca- 
tion under consideration for all the droplets in the cloud. 

If a cloud containing a block-type distribution, such as that rep- 
resented by the solid lines in figure 2, impinges upon a body the 
weighted local collection efficiency p may be defined as 


P = -i- (P^ Aw^ + ^2 ^"^2 
^t 


+ p7 Aw^) 


= P-j^ An-j^ + pg An2 + . . . + Py An-^ 


(13) 
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Equation (13) results in the 3 , value at only one surface location_ 
since 3 is a function of droplet size and surface location. If 3 
is to he determined over the entire area of droplet impingement . on the 
hody, equation (13) must he applied at enough surface points on the 
hody to evaluate 3 as a function of location s. Equation (13) is 
graphically represented in figure 4 for the surface location Sq on a 

two-dimensional hody. The value of 3 is obtained hy the summation of 
3An, which is equal to the cross-hatched area in figure 4. For other 
siirface locations (s^, Sg, S3, etc. in fig. 5) successively farther hack 

from the reference point, equation (13) is likewise applied and B at 
s^ is the area under that s^^ curve . There is a limit of impingement 

SjQ for each droplet size group (each An). This limit is determined 
hy the K and (p for each group. For surface distances, on the hody . 
greater than for any group, the 3 for that group equals zero. As 

the surface distance from the reference point increases (s > Sq) , only 

the groups of increasingly larger droplet sizes will contribute to 3> 
as shown in figure 5. By summing under each curve of figure 5, 

the ^ value for that location is determined_as previously illustrated 
in figure 4. This procedvire of summing for 3 for each location s^ 

results in a plot (see fig. 6) of F as a function of sxirface distance 
for one set of operating conditions. If a smooth droplet size distribu- 
tion c\irve is specified, equation (13) takes the form 

1 rn=l 

p = ^ / P dw = P dn (14) 

Equation (14) may be applied at several surface locations (s^) as graph- 
ically represented in figure 7, The area under any s^ curve in figure 

7 represents p for that surface location, and therefore p may be 
shorni as a function of s as in figure 8. Figures 7 and 8 for the 
smooth-type droplet size distribution curve are comparable to figures 5 
and 6, respectively, for the block-type distribution curve. 

Weighted local impingement rate Wp and weighted local collection 

efficiency p may be integrated over the surface of the body to produce 
the total impingement rate and total collection efficiency 

respectively: 

“m = f (15) 


( 16 ) 
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Equations (15) and (16) for the cloud of nonuniform droplet size corre- 
spond to equations (8) and (7), respectively, for the case of uniform 
droplet size. 


GENERAL DESCRIPTION OF EXPERIMENTAL APPROACH 

The objectives of an experimental droplet impingement study are 
the evaluations of total collection efficiency of a body, the local 
collection efficiency around the surface, and rearward extent of drop- 
let impingement on the body. The method described herein to determine 
these quantities experimentally consists of a measurement of the local 
impingement rate Wp on a body. . The measurement is obtained from a 

dye deposit left on the body. The amount of dye deposit is directly 
proportional to the water impingement rate. In this technique, water 
treated with small quantities of water-soluble dye is injected in the 
form of droplets into the air stream ahead of the body by means of spray 
nozzles. The surface of the body is covered with an absorbent material 
upon which the dyed water impinges and is absorbed. At the point of 
impact and droplet absorption, a permanent dye deposit (hereinafter 
designated as dye trace) is obtained. The amount of dye trace obtained 
in a measured time interval can be determined by colorimetric analysis 
and converted into the quantity of water which produced it if the com- 
position of the treated water is known. 

' The impingement - limit is obtained directly from the rear- 

most dye trace (where Mp-^ O). The dye analysis consists of removing 

the dye impregnated absorbent material from the body and punching out 
small segments of the material for the determination of local impinge- 
ment characteristics. The dye is dissolved out of each segment in a 
known quantity of water. The weight of dye in this solution is deter- 
mined by the amount of light of a suitable wave length that will be 
transmitted through the solution by use of a calibrated colorimeter. 

From a knowledge of the original concentration of the dye in the water 
droplets, the. weight of dye collected per unit area can be converted 
into, the weight of water which impinges at any surface location per unit 
time. This conversion assumes that no evaporation of the dyed droplets 
occurred in the air stream. The use of the experimental local impinge- 
ment rates Wp to eveiluate. Ej^ -was discussed in the previous section. 

The determination of the droplet size distribution from the experi- 
mental local Impingement rates on a body requires a thorough knowledge 
of the theoretical droplet trajectories around the body (discussed in 
a subsequent section) . A survey of the available literature indicated 
that- the most, reliable and abundant theoretical impingement data existed 
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for two-dimensional todies such as airfoils, cylinders, and rihhons. Of 
these data the cylinder analysis (refs. 1 and Z) appeared to he the most 
applicable for an experimental study that required an accurate determi- 
nation of local impingement rates from the dye traces. The impingement 
data reported herein were therefore obtained by analyzing the impinge- 
ment records from several different sized cylinders. 


APPARATUS 

Models and related equipment . - The experimental studies of droplet 
impingement characteristics on cylinders were conducted in the 6- by 
9 -foot test section of the NACA icing tunnel. Three wooden cylinders 
all of l-foot span and 2, 4, and 6 inches in diameter were Installed on 
a hydraulically actuated elevator shaft. End plates were used on all 
cylinders to maintain two-dimensional flow characteristics over the 
model. The cylinders were moved in and out of the tunnel by means of 
the elevator. A photograph of the 4-inch cylinder mounted in the center 
of the tunnel test section, is shown in figure 9.. . ' 

Pressure distributions for the 2- and 6-inch cylinders were ob- 
tained by. means of pressure belts cemented tp the cylinder sxirface.- 

A square grid made up of l/S- by l/2- by. 24-inch bar stock held 
ribbons which were used for measuring cloud uniformity and was mounted 
in the tunnel test section as shown in. figure 10. An aspirating device 
was used to determine the cloud liquid-water content and is described 
in appendix B. 

Absorbent material. - The absorbent material selected for the 
studies had to satisfy three basic requirements: (l) The material must 

have the capacity to absorb and hold at the point of impact an amount 
of dyed water sufficient to produce adequate color contrast for good 
colorimetric determinations, (2) the thickness of the material must be 
such that it can be readily formed without breaking or otherwise dis- 
torting the model contours and be of negligible thickness as compared 
with the dimensions of the model, and (3) the material must have no 
perceptible chemical reaction with the dye solution. For the studies 
herein, commercial varieties of white blotter paper of two thicknesses, 
0.017 inch and 0.023 inch, were used; however, for the models and spray 
conditions investigated no difference in results was obtained with the 
two thicknesses. . - . 

Ribbons of blotter paper l/S by 1^ inch were fastened to the lead- 
ing edge of the grid bar stock for use in measuring the cloud unl'formity 
at various locations throughout the cross-sectional area determined by 
the grid dimensions. An accurately machined paper punch was used to 
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obtain these ribbons (l/8 by 1^ in.) from sheets of blotter paper. 

Blotters about 3 inches vide were fastened to the three cylinder models 
by means of masking tape mounted parallel to the cylinder axis and dovn- 
stream of the maximum thickness point. 

Dye and water for spray fluid . - A dye that would provide impinge- 
ment records suitable for analysis was selected to satisfy the follow- 
ing criteria: (l) The dye must be readily and completely soluble out 

of the absorbent material; (2) the dye must be nonvolatile to prevent 
dye losses from the trace to the surrounding air by evaporation after 
exposure to the spray; (3) the dye must be reasonably color fast and 
chemically stable so that the trace will not fade or otherwise change 
in chemical nature during processing; (4) the solution of dye in the 
spray fluid must be a true solution and must provide sufficient color 
contrast on the absorbent material to allow accurate colorimetric analy- 
sis and determination of the limit of impingement from the rearmost 
dye trace; (5) the dye should be substantially p^ure to eliminate the 
possibility of impurities settling out of solution, and should be 
chemically inert to eliminate any reaction with the spray system 
piping. 

Of the dyes that would satisfy these criteria, a red dye of the 
azo group (Carmoisine BA Extra Cone. 118 percent factory strength. Gen- 
eral Dyestuff Corp.) was chosen for the experimental studies reported 
herein. Water-dye solutions between 0.5 and 1.0 percent dye by weight 
were studied and proved satisfactory with respect to solubility and 
color contrast. The 0.5 percent (nominal) dye in water solution was 
used throughout this investigation. A standard colorimeter was used to 
determine the light transmittal characteristics of the water-dye solu- 
tions. For colorimetric analysis with a nearly monochromatic light 
source, the best color definition for the dye-water solution was ob- 
tained at a wavelength of 520 millimicrons. 

Distilled water was' used throughout this investigation. In order 
to obtain a true solution of dye in the spray water, the dye-water solu- 
tion was agitated mechanically and heated for several hours near the 
boiling point. The addition of carbon dioxide to the distilled water 
to make the water slightly acidic seemed to improve the stability of 
dye-water solutions made with the water. 

• Spray system . - For the purposes of the present study a special 
spray system was used containing four air-water atomizing nozzle assem- 
blies of the type shown in figure Tl. The nozzle assembly ' consists of 
(l) a two-way, fast 'action, normally-closed solenoid valve, '(’2) the 
spray 'nozzle (inconel), and (3) the nozzle housing which contains the 
atomizing air inlet. These nozzles are of the same design as those used 
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in the main icing spray system of the icing tunnel and are located in 
the same low-air-velocity section of the tunnel. Each of the nozzle 
assemblies was adjustable as to its orientation, facilitating the devel- 
opment in the tunnel test section of a spray cloud that was uniform in 
liquid-water content and of sufficient size to cover completely the 
largest model studied. The dye solution was pressurized (lOO and 120 
Ib/sq in.) in a 15 gallon Inconel supply tank. From the tank, the solu- 
tion passes through stainless steel tubing and fittings to a stainless 
steel manifold and thence through rubber hoses to each individual nozzle 
assembly, (inconel, stainless steel, and rubber were unaffected by the 
dye solution; whereas other materials, including brass, copper, and 
aluminum, reacted chemically with the dye solution. ) The air used for 
the spray nozzles is pressurized (60-80 Ib/sq in.) and is nearly satu- 
rated with water vapor. This air, after passing through a porous filter 
for cleaning, is ducted to a steel manifold and then through rubber 
hoses to each nozzle assembly to be discharged through the nozzle ori- 
fice (fig. 11) at sonic velocity. 

Calibration of the spray system and valve -triggering mechanism in- 
dicated no measurable time lag in the opening and closing of the valve. 
Increasing the air-water pressure ratio from 0.5 (60 Ib/sq in. air, 120 
Ib/sq in., water) to .0.6 (60 Ib/sq in. air, 100 Ib/sq in. water) and then 
to 0.8 (80 Ib/sq in. air, 100 Ib/sq in. water) resulted in a decreased 
water flow rate. A low water flow rate for the air and water pressures 
used is accompanied by a small droplet size, whereas a large water flow 
rate is accompanied by a large droplet size. 

Air and . spray-fluid pressure were remotely and independently con- 
trolled from the main icing tunnel control room by means of suitable 
pressure regulators and transmitters. These pressures were recorded by 
means of a manometer board. All pressures were preset and could be 
repeated and maintained within ±l/2 pound per square inch. The solenoid 
valves which controlled the spray duration were linked to an electronic 
timer by means of which the duration of the spray was repeatable to 
±0.02 second. 


PROCEEURE 

Conditioning of tunnel air stream . - The studies reported herein 
were conducted at a nominal free-stream velocity of 175 miles per hour, 
a total air temperature ranging from 45° to 60° F, and a pressure alti- 
tude of approximately 28.2 inches of mercury. The tunnel air stream was 
conditioned to be as nearly saturated as possible before the body was 
exposed to the dyed-water spray. This procedure minimized the evapora- 
tion of water from the droplets during their time of travel (less than 
1 sec) from the spray nozzles to the body. (The absence of evaporation 
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was verified by weighing a ribbon on an analytical balance immediately 
after exposure and comparing the actual amount of water caught with the 
amount calculated by the dye technique.) A nearly satiorated air condi- 
tion was achieved through the control of tunnel air temperature and the 
addition of steam to the air stream. Indication of a nearly saturated 
air stream was obtained visually from the formation of a light conden- 
sation cloud in the tunnel test section prior to exposing the models to 
the dyed-water spray. Previous studies had shown that the presence of 
a light condensation cloud in the test section was accompanied by a 
relative humidity above 90 percent at the spray system; this humidity 
evaluation was made with a psychrometer . 

Exposing of models to spray cloud . - The procedure for each run 
was to preload the air and dyed-water pressure in the spray system, pre- 
set the exposure time, and then, with the tunnel air properly conditioned 
as to speed, temperature, and relative humidity, lower the cylinder with 
blotter attached into the center of the test section. In the case of 
the ribbons, the grid was mounted in the center of the tunnel during the 
air stream conditioning process. The solenoid valves in the spray nozzle 
assembly (see fig. ll) were then energized to permit the flow of dyed 
water into the air stream for the preset and automatically recorded time 
interval. The models were then removed from the tunnel. 

The exposure time for the ribbons and cylinders was varied with 
both the model size and water content in the cloud. A long exposure in 
the spray cloud (soggy blotter) causes lateral dispersion of dyed water 
in the blotters resulting in erroneous impingement rates. Excessively, 
long exposures also permit the dyed water to soak through the blotter 
paper onto the body, resulting in a loss or redistribution of dyed water. 
In the studies described herein, a satisfactory exposure- time for the 
ribbons was less than 3 seconds and for the cylinders ranged from 4 to 
10 seconds at air-water pressure ratios of 0.5 to 0.8, respectively. 

Spray cloud calibration . - Analysis of the dye collected by the 
ribbons on the grid determined the uniformity and reproducibility of the 
spray cloud. The amount of dye collected on the grid ribbons was rela- 
tively constant (±5 percent) for the portion of the cloud that covered 
a cross-sectional area in the test section greater than the frontal area 
of the largest (6-in. dlam. ) cylinder. This uniform portion of the 
spray cloud minimizes errors in local impingement rates on cylinders 
that could be attributed to a nonuniform cloud ahead of the cylinder. 

The reproducibility of the uniform portion of the spray cloud from one 
grid exposure to the next was of the order of ±3 percent variation in 
dye collected (based on average of two grid exposures). Variations of 
the cloud with respect to reproducibility were attributed to fluctua- 
tions- in nozzle discharge patterns and in tunnel air flow. 
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The uniformity of. the droplet size distribution at various loca- 
tions in the cloud, was established by placing a 2-inch-diameter cylin-^ 
der in three random positions in the cloud. The local impingement rates 
Wp for the cylinder at these positions were substantially identical; 

consequently it is reasonable to assume that the local droplet size dis- 
tribution over the entire cloud was similar. 

Dissolution of dye from absorbent material . - Before dissolving the 
dye from the blotter paper for colorimetric -processing, the l/S- by 

l^inch ribbons from the grid were accurately trimmed to 3/4 inch in 

order to eliminate any end effects caused by the masking tape. The 

punch used to obtain the l/8- by li-inch ribbon blotters was used to ■ 

remove segments from the cylinder blotters for local dye -impingement 
determinations. No trimming of these cylinder blotter segments was re- 
quired because the punched segments are sufficiently inboard from any 
possible edge effects. A typical cylinder blotter with punched segments 
removed is shown in figure 12. The punched blotter segments were placed 
in individual 15 milliliter capacity test tubes and the dye dissolved 
out of each segment by adding a known amount of distilled water suffi- ■ 
cient to completely cover the blotter segment (generally 5 ml). The 
dissolution of dye from the blotters was expedited by placing the test 
tube containing the distilled water and the dyed blotter segment in a 
water bath and heating for about 10 minutes at a temperature slightly 
below the boiling point. Mechanical agitation of the test tube was not 
employed, since agitation would promote excessive disintegration of the 
blotter paper. The test tubes were each covered with a vented cork to 
prevent the loss of a significant volume of water during heating. After 
heating was completed, the dye-water solutions were allowed to cool to 
room temperature and the effective concentration of the dye solution was 
determined colbrimetrically . ' It was ascertained from successive dilu- 
tions of the same blotter segment that the initial dilution dissolved 
into solution more than 99 percent of the dye contained in the segment; 
therefore, one dilution was considered adequate. 

Colbrimetri-c analysis . - After the effective concentration of the 
solution in each test tube was determined from the colorimeter, a correc- 
tion was made because the' solution in each test tube contained a small 
amount of blotter fiber in suspension. These fibers absorbed some of 
the light in the colorimeter, causing more total light' to be absorbed 
by the contents of the cell than- would be absorbed if only a piore solu- 
tion (no- fibers)- were analyzed-. The true solution concentration P is 
the effective- concentration of the solution minus the effective concen- 
tration of a- •"clean -blotter" ' solution. A clean blotter solution is one 
which contains the- same amount of water as the other solutions, but the 
blotter segment therein has no dye on it. Since the amount of blotter 
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fiber in suspension appeared to be a function of the heating time, sev- 
eral clean blotter solutions were heated along with each batch of test 
tubes. An average of the effective concentrations of these clean blotter 
solutions was considered to be the best correction for that particular 
batch. Each batch consisted of all the solutions obtained from one cyl- 
inder blotter (19 from the blotter shown in fig. 12). 

The amount Mp of dye impinging on any blotter segment per hour 
per square foot is given by: 


Mp (lb/ (hr)(sq ft)) 


t 

36W 


Fb 

AAg 4.536x10^ 


7.94x10"^ 


Pb 

tAAs 


(17) 


The weight of water (which is approximately the weight of solution) per 
unit weight of dye is a function of the concentration c, the percent by 
weight of the original dye-water solution placed in the spray system, so 
that 


M- 100 
P 



c 


(18) 


Equations (17) and (l8) result in 


Wp (lb/(hr)(sq ft)) = (19) 

s 

The local water impingement rates given by equation (19) are the 
primary results obtained by the dye-tracer techniques and are believed 
to be accurate and reliable (see RESULTS and DISCUSSION sections). How- 
ever, the correlation of the experimental impingement characteristics 
obtained from these Wp values with the independent parameters (K, <p, 

etc.) requires an accurate knowledge of the droplet size distribution. 
Because of the inadequacies of presently available instruments for de- 
termining droplet size distribution, a method has been developed by 
which the droplet size distribution can be ascertained. This method, 
which is presented in the following, section, permits the calculation of 
the droplet size distribution by relating the experimentally obtained 
Wp values, for a body to the theoretical trajectory results for the 

same body. 


METHOD FOR DETERMINING DROPLET SIZE DISTRIBUTION 

The generalized dependent impingement parameters discussed in the 
section entitled BASIC IMPINGEMENT RELATIONS for both uniform droplet 
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size .and nonuniform droplet size distributions are simplified for the 
cylinders used herein because the cylinders are symmetrical and are two- 
dimensional bodies. The resultant equations for a cylinder are: 



( 20 ) 


= 0.658 Vo EE^ 


( 21 ) 


P = lim 
As-K) 


Ayo <iyo 


As 


ds 


( 22 ) 


The surface distance s on a cylinder. is measured from the stagnation 
point and may also be denoted in terms of the central angle 9 measured 
from the stagnation point, while the limit of impingement is denoted by 


The term P may be determined from the theoretical trajectory re- 
sults for any given set of impingement ccJnditions (K, <p) by plotting 
y^/R (trajectory starting ordinate at infinity) as a function of s/r 

(angle of impingement, radians), as shown in figure 13 (obtained from 
ref. l). The slope of this cvirve (see eq. (22)) at any value of s/r 
yields the value of P at that point. A 'series of curves of P (ref. 
l) are shown in figure 14 plotted as a function of 9 for several sets 
of impingement conditions (combinations of K and <p). An average, of 
the <p = 1000 group and the <p = 10,000 group has been cross-plotted by 
the use of figure 15 (from ref. 1^) resulting in , (where p is 

expressed as a function of 9 and Qj^) as shown in figure 16. The 

curves in figure 16 are essentially independent of <P for 1,000 ^ <p 
< 10,000, and are required later in the analysis of the experimental 
impingement data.- 


^In obtaining figure 15 of the present report from figure 7 of 
reference 1, the cxjrves of figure 7 were extrapolated to pass through 
0JJI = 0 at K = l/S. In this extrapolation, the existence of a critical 

value of K (K^. = l/8 for cylinders) was assumed. Reference 2 

states that, for conditions specified by values of K less than K^, 

no droplets, even those approaching along the air stagnation line, will 
impinge on the body; and if all impingement conditions except droplet 
size are specified, there is a critical 'droplet size af. which, will 
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Langmuir suggested (ref. 13, p. 4) that it might he possible to 
determine droplet size by comparing theoretical trajectory results for 
cylinders with experimental measurements of rime ice thickness on a 
cylinder. If theoretical trajectory results (such as those in refs. 1 
and 2) were compared with values obtained from measurements of ice 

thickness, the results would be of doubtful value, since the presence 
of the ice would change the effective contour of the cylinder and in- 
validate the applicability of the theoretical results. The dye-tracer 
method does not involve any appreciable change in cylinder contovir and 
consequently is a more reliable sovirce for Wg values. In the follow- 
ing method, the experimental Wp values for a cylinder are combined 

with the theoretical curves of figures 15 and 16 to determine the drop 
size distribution for a given set of impingement conditions (where <p 
is known and K is a known function of the droplet size). The basic 
equation which must be solved to determine the droplet size distribu- 
tion is obtained by multiplying both sides of equation (14) by 
0.329 Vq w^ which, with equation (lO), results in 

/ w=w^ 

( 23 ) 

The subscript bracket notation indicates that Wp is a function of 9 
and that P is a function of both 9 and (fig. 16). Each value 

of droplet radius determines a value of K for a cylinder which^ with 
specifies a value of a limit of impingement 9 ^ (fig. 15); hence^ 

the problem of obtaining droplet size as a function of water content w 
is reduced to the problem of finding w as a function of 0^. The 


result in K = K^. It is shown in references 2 and 3 that varies 

with body shape^ depending primarily on the flow field near the stagna- 
tion region. Because of the physical dimensions of the droplet (the 
droplet is considered to be a point in the theoretical analyses)^ the 
actual existence of this critical condition appears to be questionable^ 
and no rigorous theoretical or experimental proof of the physical phe- 
nomenon occurring when K < has been presented. The approximate 

equation of motion of the droplet does change form^ however^ when K = K^. 
When the accuracy of the theoretical data for small values of K is con- 
sidered, it appears reasonable to extrapolate the curves to agree with 
the critical K concept. For the analysis presented herein, it has been 
assumed for simplicity that there is a critical K equal to l/S for cyl- 
inders. This assumption does not alter the results appreciably from those 
which would be obtained by ass\iming that ®m approaches zero asymptoti- 
cally as K approaches zero in figure 15. 
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exact solution of equation (23) for 0^ as a function of w is a dif- 
ficult problem in the calculus of integral equations because the unknown 
function is neither nor Pj"g It is necessary to determine 

a relation between 9.^ and w such that equation (23) will be satis- 
fied for any value of 9. A method has been developed by which a close 
approximation to the exact mathematical solution may be obtained. The 
following steps describe this method, and the method is paralleled by 
a nttmerical example in appendix C: 

Step 1: Divide a typical experimental Wp against 9 curve for 

a cylinder (solid line, fig. 17) into eirbitrary increments of 0,_as 
shown by the vertical dotted lines in figure 17. The values of Wp are 

obtained from the dye trace on segments taken from blotters such as 
shown in figure 12. Also shown in figure 17 are the contributions to 
Wp of each droplet size group (dashed curves), the determination of 
which will be discussed in the following paragraphs. 


Step 2: Assume that the unknown w against 0j^ curve is composed 

of a block-type droplet size distribution like the dashed lines in fig- 
ure 18, where each droplet size group has a 0^ equal to one of the 

arbitrary 0 values assigned in the first step, that is. 


0m. = 0 


“max 


= 0 . \ 




The value of Aw for each group, however, imst still be determined. 

Also shown in figure 18 is a solid curve representing a fairing of the 
block-type distribution. 

Step 3: The Wp value obtained by analytically weighting a block- 

type droplet distribution curve for w against 0^j ■ must be equal to 
the experimental Wp value at each of the chosen 0 values. 

The third step determines the amount Aw^^ of water contained in 
droplets of size a^^ specified by 0jjj^ because all water hitting from 
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02 ^ (fig* 17) is made up of droplets of size specified 137 0,^^ 

(condition 2, preceding). Therefore, at 03 where only droplets of 

size specified l>y ^m^^ impinge, equation (23) may "be used and ^^0j2j,8j 

becomes the constant Ppg g n and w^ is Aw^: 

L mi^ 2J 

/ w=Aw2^ 
w-0 


= 0.329 Vq P[-0„^^,02j ^'"l 


or 


Awi = 




1 - 0.329 Vo 


(24) 


Thus the magnitude of Aw^ is determined and the values 0j^^ and Aw-j^ 

describe the first droplet group, which is indicated as the bottom block 
in figure 18. 

After the value of Aw-j^ has been determined, the Wp value con- 
tributed by the first group characterized by 0j^^ and aj_ may be deter- 
mined for all values of 0, since equation (23) shows that 




(25) 


where ej 1® determined from figure 16. ^p^0ui^,0j 

indicated as the height of the bottom dashed curve (curve l) in figure 17 


The water content Aw2 for the next droplet size group now may be, 
determined, since for 0j ^ 0 02 only droplets of group (l) and group 

(2) hit. At 0 = ^3 the amount of Wp due to group (l) has already 
been determined by equation (25). Therefore, at 83 the amount 
Wp^0 due to group (2) is 
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Wo 


'Ppm2^03] = (26) 

as indicated in figure 17. This value of Wc 


since 


^3 *0jjj determines Aw2, 


2 0.329 Vq P [ 0^12,03] 


Awo = 


(27) 


The value of 3 again is determined from figure 16 at 0jjj 

0 = 03. The second group may then he indicated as another 
figure 18. 



block 


and 

in 


The value of Wp due to the second group may he determined for all 

values of 0 and indicated as the difference in height between dashed 
curves 1 and 2 in figure 17, since 

«?[9M2,e] = ''0 ^2 P[em2.9] 

where 3 is obtained from figure 16. 

The amount of water impinging at 0^ due to the third group may 

then be determined because the amounts of water contributed by the first 
and second groups of droplet sizes are known. 

The value of Awj is then determined by an equation similar to equa- 
tions (24) and (27). This process is continued until a value of Aw 
has been determined for each droplet size group. Each group with its 
contribution to 'the total water content is indicated as a block in 
figure 18. 

The block diagram of figure 18 is an approximation of the droplet 
size distribution (in, terms of 0jjj) to which the cylinder was exposed. 

This approximation is inaccurate to the extent that the weighted Wp 

(the sm of the heights between the dashed curves, fig. 17) equals the 
experimental Wp only at the selected 0 points. Small Increments of 

0 may be taken to cause the weighted Wp curve to approximate the 
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experimental Wp curve to any degree of accuracy desired; however, 

satisfactory results may he obtained hy a faired curve based on only 5 
or 6 increments, as shown by the solid curve in figure 18. This faired 
cxxrve is reasonably independent of the choice of 9 increments and 
when it is used to weight Wp, as described in a previous section of 

this report, a very close approximation to the experimental Wp curve 
is obtained. 


DETERMINATION OF CLOUD LIQUID-WATER CONTENT 

As previously discussed, the data from the solid curve of figure 18 
may be replotted in terms of droplet radius and cumulative liquid-water 
content as shown in figure 19. The droplet size distribution cxirve in 
figure 19 is defined only for droplet sizes larger than a^; the value 

of ac depends on the diameter of the cylinder used. Completion of the 
droplet-size distribution curve through the use of progressively smaller 
cylinders (with accompanying smaller a^ values) is indicated. This 

procedure would define a greater portion of the droplet size distribu- 
tion and also result in a value of w at which would approach the 

total liquid-water content of the cloud w^. Because it is impractical 

to obtain local water impingement rates on very small cylinders with the 
dye-tracer technique, the total liquid-water content w^ was obtained 

by a simple aspirating device. Details of this device and its use in 
obtaining the total liquid-water content are presented in appendix B. 


RESULTS 

The experimental data are presented in terms of local rates of 
water impingement Wp as a function of cylinder central angle 0. 

These data were analyzed to obtain the droplet size distribution of the 
impinging cloud by, use of the method outlined in a previous section. 

Local impingement rates and impingement limits . - The local im- 
pingement rates Wp for the three cylinders as determined from the dye 

traces by use of equation (19) are shown in figure 20 as a function of 
the cylinder central angle 0 for three air-water pressure ratio set- 
tings of the spray nozzles . It is apparent from this figure that for a 
given droplet size distribution (fixed air-water pressure ratio), the 
smallest cylinder yields the highest local impingement rates. The local 
impingement rate for each cylinder also increases with increasing drop- 
let size (as indicated by the decreasing air -water pressure ratios). 
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The limits of impingement (Wp = 0) on the cylinders for the spray 

conditions studied are also shown in figure 20. It is apparent that as 
the cylinder size is decreased, the impingement on a cylinder occurs 
over a greater cylinder central angle 0. The average limit of impinge- 
ment for three repeat runs as defined hy the central angle ®m^^ varies 


with the air -water pressure ratio. For the 6-inch cylinder the limit of 
impingement changes from 46° to 74° for a change in air-water pressure 
ratio from 0.8 to 0.5. Comparable values for the 2-inch cylinder are 
70° and 84°. 


Maximum droplet size. - The maximum droplet size is determined from 
the experimental impingement limit (Wp -► 0, fig. 20) hy the relation of 

0^, K, and <p shown in figure 15. With- (p and 0^^ known, the droplet 

size can he readily calculated from the K parameter. Although the 
TnRviwniTn droplet size obtained from different cylinder sizes at the same 
spray condition varies slightly (the smallest cylinder yields a somewhat 
larger droplet size value than the largest cylinder), droplet size dif- 
ferences are well within the range where the experimental impingement 
limit can he determined. At large values of the slopes of the 

theoretical K against 0^ curves (fig. 15) indicate large changes in 
K for, only small changes in 0jj^j hence, a small error in 0j^ will 

magnify the error in K or droplet size. In general, therefore, the 
ma-iciTmiTn droplet size should he specified by cylinders having the maxi- 
mum angles of impingement in a , range where the slope of the K against 
0JJ ctirve is steep (i.e., the maximum droplet sizes determined hy the 

largest cylinders used herein are more reliable than those determined 
hy the smallest cylinder). The maximum droplet sizes, based on refer- 
ence 1, as a function of the air-watef pressure ratio are summarized as 
follows: 


Air -water 

pressure 

ratio 

Maximum droplet radius, 
a, microns 

(average for 3 repeats 
on 6 in. cylinder) 

0.5 

28.0 

.6 

20.4 

.8 

■■ ■ 11.9 . 


An analysis using the K against curves of reference 2 indicates 

a maximum droplet size approximately 15 percent greater than that ob- 
tained with the data of reference 1; however, this percentage difference 
in droplet size is generally in accord with the differences existing 
between these two theoretical studies. 
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Droplet size distribution and total liquid-vater content . - The 
cumulative liquid-water content contained in the spray cloud may he de- 
termined as a function of droplet size from the results of figure 20 by 
the procedures discussed in the section METHOD FOR DETERMINING DROPLET 
SIZE DISTRIBUTION. A detailed example of the numerical procedures is 
given in appendix ^ The curves of cumulative liquid-water content w 
as a function of .droplet radius a based on the theoretical curves of 
figures 15 and 16 (ref. l) for all three cylinders are shown in figure 
21 for air-water pressure ratios of 0.5, 0.6, and 0.8. In figure 21, 
in the large droplet size range, the scatter among the results from the 
three different cylinder sizes at a constant spray condition is of the 
same order as the scatter between repeat runs of the same cylinders. 

For small droplet sizes the separation of the results relative to cyl- 
inder size (hereinafter referred to as ’’cylinder-size trend”) is quite 
evident. This cylinder-size trend should not exist since, at constant 
spray conditions, the cloud has one droplet size distribution. Similar 
w against a curves axe shown in figure 22 based on the theoretical 
data of reference 2. In general, these' curves (fig.' 22) are in agree- 
ment with those obtained using. reference 1 (fig. 21); however, the drop- 
let sizes are slightly larger and the cylinder-size trend in the small 
droplet size region is not so apparent. The differences noted (between 
figs. 21 and 22) are again caused by the differences inherent in the 
theoretical results of references 1 and 2. Because more data points 
are presented and the accuracy of the data is stated, reference 1 will 
be used hereinafter rather than reference 2 as a' basis for the subse- 
quent results. 

In obtaining a single curve which is the best representation of 
the actual droplet size distfibiition existing in the tunnel for a given 
set of spray conditions, it is reasonable to use an average of the nine 
curves (three repeats times three cylinder sizes) in the larger droplet 
size range where the results from all three cylinder sizes are in 
reasonable agreement. In the range of the smaller droplet sizes (a < 6 
microns, see fig. 2l), it appears that the theoretical data are not 
strictly applicable because a different droplet size distribution curve 
is obtained for each cylinder size. However, the end point of the drop- 
let size distribution curve (a = 0, v = v^) is single- valued and is 

specified by the values obtained with the aspirating device. The 

averages (for several repeats) of the total liquid-water content w.^ 
obtained with the aspirator are summarized as follows: 
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Air-water 

Liquid-water content. 

pressure 

'^t 

ratio 

g/cu m 

0.5 

0.70 

.6 

.58 

.8 

.43 


Because of the cylinder-size trend for droplet radii less than approxi- 
mately 6 microns (fig. 21), a straight line was faired between a *6 
microns and the tabulated value of w^ for each pressure ratio. The 

results of this interpretation of the data are shown in figure 23; the 
reasoning underlying this interpretation is presented in the DISCUSSION 
section. 


Volume median droplet size and dimensionless distribution curves. - 
One of the more useful impingement parameters in correlating droplet 
size data is the volume median droplet size (ref. 2). For the three 
pressure ratios given, the following volvime-median droplet sizes are ob- 
tained from figure 23: 


Air -water 

Volume median droplet 

pressure 

radius^ a, microns 

ratio 


0.5 

7.4 

.6 

6.0 

.8 

3.8 


Droplet size distributions are often plotted in terms of the dimen- 
sionless ratio of droplet size to the volume -median droplet size 

as a function of percent cumulative volume of water in the cloud. A 
plot of a/ajjjg^ as a function of the percent cumulative volume n = w/w^ 

for the three air-water pressure ratios is shown in figure 24. Although 
the experimental curves are not identical, the deviations are believed 
small enough to warrant the approximation that the three air-water pres- 
sure ratios used produce substantially the seuae dimensionless droplet 
size distribution. 


DISCUSSION 

The impingement data obtained by use of the dye-tracer technique 
are considered valid; however, certain discrepancies which cause devia- 
tions between the exp>erimental results and theoretical analyses require 
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clarification. The apparent cylinder-size trend, shown in figure 21, 
should not occur for a given spray condition since all the cylinders are 
exposed to the same cloud. Probable reasons for the cylinder-size trend 
appear to be a combination of (l) the effect on droplet drag of droplet 
acceleration (ref. 14) and (2) the difference between the flow field 
about a cylinder assumed in references 1 and 2 (ideal fluid flow) and 
the flow field existing in the experimental studies. (Refs. 1 and 2 
assume that the drag coefficient is a function of only the relative 
Reynolds number; however, reference 14 shows that the drag coefficient 
is also a function of the rate of change of the relative Reynolds number. ) 
The effect on droplet drag of droplet acceleration in the flow field is 
to increase the drag coefficient, resulting in reduced local impingement 
rates, total collection efficiencies, and limits of impingement. This 
effect will vary with body size because a droplet experiences the accel- 
eration for a longer time while in the flow field caused by a large, body; 
hence, the effects may be greater for large bodies than small bodies. 

The effect of acceleration drag should also be less noticeable for large 
droplets which are not accelerated so much as small droplets. Conse- 
quently, acceleration drag would tend to produce a cylinder-size trend 
shown in figure 21. 

A comparison of the theoretical and experimental surface velocities 
on cylinders is shown in figure 25. It is apparent that the local vel- 
ocity ratios differ over the entire cylinder surface. The early separa- 
tion of the flow (0 = approximately 50°) and the deviation of velocity 
ratio from the theoretical values (0 < 50°) are typical of bluff bodies 
studied in other wind tunnels (ref. 15). The difference between experi- 
ment and theory is due to the presence of a bluff body in a nonideal 
fluid and is possibly magnified by the tunnel air-stream turbulence 
level. These considerations indicate that such large differences in 
local velocity on the cylinder may mean that the flow field ahead of the 
cylinder is also considerably different from the theoretical flow field. 
Consequently, the actual droplet trajectories may be different from the 
theoretical calculations of reference 1 and may contribute to the 
cylinder-size trend shown in figure 21. In view of the possible effects 
of acceleration drag and flow field deviations, the agreement of the ex- 
perimental data shown in figure 21 appears reasonable. 

Additional factors that have been considered to account for the 
cylinder-size trend include: droplet evaporation effects, . splatter of 

droplets on blotter paper, effect of blotter texture on local water 
catch, oscillation or pulsation of spray cloud in the air stream, initial 
relative velocity differences between droplet and air, droplet distortion 
effects on droplet drag coefficient, boundary layer considerations, and 
accuracy (with regard to the mechanics of the analyzer) of the theoret- 
ical trajectories (ref. l). Evaluations of the effect of these factors 
on the cylinder-size trend all indicated a negligible effect or an effect 
■that would produce a minor trend opposite to . that obtained in the experi- 
mental study. 
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As stated in the RESULTS and presented in figure 23, a straight 
line was assumed to define the droplet size distribution curves between 
droplet radii of approximately 6 and 0 microns. The extremities of this 
fairing appear reasonably well defined: (l) Analysis of droplet size 

values from all cylinders in figure 21 are in general a,greement for 
a > 6 microns and (2) at a equal to zero in figure 23, the value of 
total liquid-water content from the aspirator is considered valid. Be- 
tween droplet radii of 6 and 0 microns (see fig. 23), three general fair- 
ings are possible: (l) a straight line as shown, (2) a concave fairing, 

or (3) a convex fairing. The concave fairing would indicate that the 
cloud contains a multitude of extremely small (of the order of a = 1 
micron) droplets and very few droplets of the or^r of radius 2 to 5 
microns. The convex fairing, on the other hand, would indicate a large 
number of droplets of the order of 2 to 5 microns in radius and very few 
droplets less than about 2 microns in radius. 

As previously shown in figure 21, the droplet size distribution 
curve (obtained with the theoretical trajectory data of ref. l) ap- 
proaches the concave-type fairing for droplet sizes <6 microns in 
radius; however, each cylinder size indicates a considerably different 
concave fairing of droplet distribution data for the same spray condi- 
tion. This contradiction makes any fairing based on these theoretical 
trajectories unreliable in this droplet size range. 

It should be possible to obtain values of total liquid-water con- 
tent from the experimentally detemiined lnq>ingement rates on the grid 
ribbons by using the theoretical trajectory data of reference 2. If the 
w^ values so determined are to agree with the w^ values obtained with 

the aspirating device, the droplet size distribution curve (in the range 
of droplet sizes <6 microns in radius) would necessarily require a con- 
cave fairing. However, the theoretical trajectory results for a ribbon 
(ref. 2) are subject to the same limitations (such as neglecting accele- 
ration drag etc.) which are inherent in the cylinder trajectories (ref. 

1) and which are believed responsible for the cylinder-size trend shown 
in figure 21. The theoretical trajectory data for ribbons and cylinders 
therefore are not considered an adequate guide in fairing the droplet 
size distribution curves of figure 23 for droplet radii less than approx- 
imately 6 microns. 

In the small droplet size range (a < 6 microns), unpublished data 
from oiled slides and a droplet sampler currently being developed, al- 
though subject to sampling errors due to unreliable values of collection 
efficiency, show no Indication of a concave-type droplet size distribu- 
tion. These data indicate a convex shape or possibly a straight-line 
fairing. 

In general, an abrupt lack of droplet sizes of the order of 2 to 5 
microns in radius coupled with a multitude of droplets less than 1 micron 
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in radius appears to te iroprotatie . On the other hand^ a distribution 
of droplet sizes from 6 microns in radius down to some finite hut limit- 
ing size such as 1 or 2 microns in radius (convex distribution) appears 
feasible with an air-atomizing nozzle of the type used herein. 

Because all available evidence supporting the concave-type fairing 
for the range of droplet radii less than approximately 6 microns in fig- 
ure 23 appears contradictory and because the experimental evidence sup- 
porting the convex fairing is insufficient to determine the extent to 
which the curve should be convex, the straight-line (dashed) fairing of 
figure 23 is used. 

The accuracy of determining droplet size distribution by the use of 
the dye-tracer technique together with theoretical trajectory data could 
be definitely established if an accurate instrument were available to 
measure directly a representative sample of all droplet sizes present in 
the experimental spray cloud. The present state of development of the 
camera (ref. 16) or the oiled-slide technique (ref. 17) for measuring 
droplet size is not considered so adequate a means for determining drop- 
let size as the indirect method herein described. 


CONCLUDING REMARKS 

A technique for experimentally determining the impingement charac- 
teristics of an arbitrary body and a method whereby the droplet size 
spectrum of the spray cloud can be determined from the experimental im- 
pingement on a body for which the theoretical trajectories are known 
have been developed. These developments present an opportunity to study 
experimentally the impingement characteristics of other shapes such as 
airfoils, radomes, inlet cowlings, and asymmetrical bodies. Although 
analytical trajectory studies can be made for these various bodies, the 
problems associated with the determination of the flow fields and con- 
sequent droplet trajectories often prove difficult and tedious. Use of 
the dye-tracer technique for arbitrary shapes in any suitable testing 
facility would proceed as described herein for cylinders: The droplet 

size distribution in the spray cloud would be established using bodies 
for which the theoretical impingement characteristics are known, while 
an aspirating device similar to the one described (appendix B) may be 
used for the determination of total liquid-water content. The unknown 
body with blotter then would be exposed and analyzed for local rates of 
water impingement, including the determination of the maximum extent or 
limit of impingement, and integrated for total collection efficiency. 

Impingement limits 'and collection efficiency are often presented 
in terms of K and q> (where K is indicative of the inertia of the 
droplet and (p represents the deviation of the drag forces from Stokes' 
law) for correlating impingement characteristics. Langmuir has suggested 
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the use of a modified K parameter. Reference 18 discusses and illu- 
trates previous analytically determined impingement data (for cylinders, 
airfoils, ribbons, spheres) in terms of a Kq parameter (Kq = (X/Xg) K) 

where X is the true range a droplet would have when projected into 
still air and Xg is the range of a c^oplet obeying Stokes law. Plot- 
ting a dependent impingement parameter, such as total collection effi- 
ciency, as a function of Kq yields an essentially single -valued curve 

(for most of the simple bodies for which trajectory results are presently 
available) independent of <p, this curve being the Stokes’ law solution. 
The extension or interpolation of experimental as well as analytical data 
points over a wide range of the pertinent impingement variables - droplet 
size, body size, speed, and altitude - is greatly facilitated by this Kq 

parameter, even though no rigorous theoretical proof of the significance 
of Kq is available at this time. 

The dimensionless impingement parameters previously discussed (K, 

<p, \lf, RCq, and Kq) are presented in the literature in terms of clouds 

containing droplets of the same size. Analytical calculations have shown 
that the presence of a droplet size distribution does not alter the use- 
fulness of the Kq parameter if it is evaluated with a pertinent cloud 

characteristic. Using volume-median or -maximum droplet size of the dis- 
tribution to calculate Kq will tend to reduce total collection effi- 
ciency or limit of impingement, respectively, to single-valued curves 
(Stokes' law solutions). The total collection efficiency curve as a 
function of (KQ)jjjg^j^ for a droplet size distribution will be somewhat 

skewed from that obtained with a uniform droplet size spectrum; the 
amount of skew depends on the droplet size distribution. Reasonable ex- 
tension and interpolation of experimental data for conditions other than 
those studied is possible with the Kq parameter. 

.Further analysis has shown that the dependent impingement parameters 
for the case of uniform droplet size (such as obtained from trajectory 
results )_may be obtained from a series of experimental impingement rate 
curves (Wp against s/l) when the experimental droplet size distribution 

is known. The method requires curves of Wp against s/L, which may be 

obtained by the dye-tracer technique for either (l) one body in the pres- 
ence of several (say about six) different, but known, droplet distribu- 
tions or (2) several sizes of the body in the presence of one droplet 
distribution. The method of calculation is a step-wise approximation 
and may be derived by methods similar to those used for obtaining the 
droplet distribution in the section METHOD OF DETERMINING DROPLET-SIZE 
DISTRIBUTION.. The results of, such an analysis would be useful in (l) 
calculating the impingement rates for a body in a cloud of. arbitrary 
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droplet size distribution, and (2) determining whether the assumptions 
used in theoretical trajectory calculations are valid. This determina- 
tion requires an independent means (reliable instrument) for measuring 
the droplet size distribution. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 3, 1954 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 
area, sq ft 

area (perpendicular to air stream) at x = - sq ft 
droplet radius, microns (3. 281x10“® ft) 

volume of water used to dissolve dye from blotter segments, ml 

percent concentration by weight of dye in water solution used in 

spray system — x lOO percent s ^ 100 percent 

lb solution lb water 

total collection efficiency in clouds of uniform droplet size 
defined by eq. (1), dimensionless 

2.19x10“^® a^VQ 

inertia parameter, , dimensionless 

modified inertia parameter, K, dimensionless , 

^s 

a significant body dimension: chord length of airfoil, radius of 

cylinder, one-half ribbon width, ft 

local dye collection rate, lb/(hr)(sq ft) 

cumulative fraction of total liquid-water content = w/w^, 
dimensionless 

true concentration of solution obtained from blotter segments, 
mg dye 
ml solution 

cylinder radius, ft 
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R6q free-stream Reynolds number with respect to droplet, 

9.62x10-6 apVo , 

= YK <p, dimensionless 

s distance along sxxrface from point of reference, ft 
t e3q>osure. time, sec 

Vq free stream velocity, mi/hr 


W_ total water impingement rate in cloud of uniform droplet size, 
lb/(hr) 

Wo local water impingement rate in cloud of uniform droplet size, 
lb/ (hr) (sq ft) 

w cumulative liquid-water content contained in droplets of sizes 
from ajnax particular droplet size, g/cu m 

w^ total liquid-water content of cloud, g/cu m 


X abscissa of coordinate system parallel to free-stream velocity 
direction, ft (see fig. l) 

Yq ordinate of coordinate system perpendicular to free-stream vel- 
ocity direction at x = - <», ft (see fig. l) 

P local collection efficiency in cloud of uniform droplet size, 
defined by eq. (4), dimensionless 

9 central angle of cylinder measiired from x-axis, 57.3 s/r, deg 

X true range of droplet as a projectile injected into still air, ft 

\q range of droplet as a projectile following Stokes' law, ft 

(i viscosity of air, lb/(ft)(sec) 

p density of air, Ib/cu ft 

0.423 p2VqL 

<p Independent impingement parameter> , dimensionless 
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ilf Independent impingement parameter^ ^: . ^9 dimensionless 

Subscripts: 

c critical condition for impingement 

f frontal, projected parallel to free-stream velocity direction 

m tangent trajectory 

med volume median 

max maximum droplet size 

s surface 

particular values of a variable 

2,3,n 

^ ] brackets denote functional notation 

Superscript: 

— weighted value due to effects of more than one droplet size 
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APPENDIX B 


DESCRIPTION OF ASPIRATING DEVICE 

An aspirating device was used to obtain measurements of the total 
liquid-water content of the cloud. This device consisted essentially 
of a tube which sucked in the approaching air and cloud droplets at the 
free -stream velocity (inlet velocity ratio^ l.O) so that both the air 
streamlines and droplets entered the tube along straight-line paths. 

All the dyed droplets were deposited on a filter mounted within the 
tube^ leaving a dye trace which could be analyzed color imetrically. 

The inlet tube consisted of stainless steel tubing (10- by l/4-in. 

0. D. by 0.014-in. wall thickness) with the front end tapered to the inter 
nal diameter. Behind the inlet tube was a diffuser section increasing in 

1. D. to 3/4 inch. The increase in cross-sectional area decreased the 
velocity through the diffuser, and consequently reduced the pressure 
drop in the following section where the filters were mounted. The 
filters consisted of small circular cloth pads spaced about I/I 6 inch 
apart axially w A sufficient number of filters were used so that after 
exposure the rearmost filter was always free of dye, indicating that no 
dye passed downstream of the filter section. Air flow through the aspi- 
rating device was obtained by the use of an ejector attached aft of the 
filter section with an exhaust line. .The suction in this exhaust line 
could be regulated to adjust the weight flow through the tube. The 
pressure was measured with a static tap on the inside of the inlet tube 
3/4 inch from the front end. The tap had been calibrated to give pres- 
sure as a function of the weight flow through the tube and was used as 
an indicator for obtaining a weight flow equivalent to an inlet velocity 
ratio of 1.0. The device was mounted in the tunnel with the inlet tube 
parallel to and facing into 'the air stream. The suction in the exhaust 
line was adjusted to obtain the correct inlet velocity ratio, and then 
the device was exposed to the spray cloud in the same manner as for the 
cylinder models. After exposure, the filters were removed, and the dye 
in the filters along with any dye left on the inside of the inlet tube 
and diffuser was dissolved into a known quantity of water. The concen- 
tration of the resulting solution was measured colorimetrically. 

Equation (19) was modified to evaluate the rate of water catch for 
the aspirating device as follows: 

^ 0.794 

= 5— 

tc AAg 


where 
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Wp total rate of water catch W^jj divided hy frontal area Af of in- 
side cross section of tube 

P* effective concentration (mg/ml); also considered true concentration 
because only a negligible amount of filter fibers were in the 
solution 

b* total amount of water (ml) used to dissolve dye from inside walls 
of device and filters 

AAg eqml to frontal area Ap of inside cross section of tube 


Because all the droplets directly ahead of the inside cross- 
sectional area of the tube £nter the tube (straight-line trajectories), 
the collection efficiency Ej^ of the aspirating, device is 1.0, and the 

Wp value may be used in equation (9) to determine the total liquid- 
water content of the cloud: 


^t 


0.329 VoAfE^ 


0.329 VoE^ 


0.794 P*b* 
0.329 tc AAgVo 1.0 
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APPENDIX C 


NUMERICAL EXAMPLE OF DETERMINATION OF DROPLET DISTRIBUTION FROM 
EXPERIMENTAL LOCAL WATER IMPINGEMENT RATE 

The Wp values for this example are obtained from the repeat curve 

for the 2-inch cylinder in figure 20(a). For these data, K = 0.0387a^ 
and (p = 2800. 

A convenient setup for use in calculating the droplet size distri- 
bution is shown in table I. The first column represents the arbitrary 
values of 8 at which the Wp curve was divided. The second column 
represents the corresponding values of Wp (i.e., Wpj-g^j = 2.99). The 

top row represents the values of each equal to a 8 value and each 

specifying a droplet size group. In any box, determined by the inter- 
section of a 0 row and a 0^ coluunn, the number above the diagonal 

line is the value P for that value of 0^ at that value of 0 as de- 
termined from figure 16 (i.e., P^ 0 m^, 07 j = P[40°,20°] “ 0.335). The num- 
ber below the diagonal line in any box is the corresponding value of Wp 
(which is part of Wp at that value of 0). The object of the calcula- 
tion is to determine the values of Aw in the second row. The procedure 
outlined in the METHOD OF ANALYSIS section gives the following results: 

Group 1: • 

"P[V®2] ' “PPa] ' 


■ 0.329 V(j e 


"P[%.82] 
® ''O 


0.62 

57.6(0.237' 




Beneath diagonal line 
in 0^^ column 


= (57.6)(0.0454)(0.237) = 0.62 for 0p 
= (57.6) (0.0454) (0.406) = 1.06 for 63 

= ( 57 . 6 ) (0.0454) (0.942) = 2.46 for 0g (25) 
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For subsequent droplet size groups the numerator of the equation for 
AWj^ is the Wp at minus the sum of the values contributed 

by the preceding droplet size groups in the particular 0^+1 

Group 2: 


1.93 

^“2 " 0.329 Vg ' (57.6)(0.183) “ 



= 57.6xAw2 



Beneath diagonal line 
in 0™ column 


= 57.6 (0.183) (0.183) = 1.93 for 63 

= 57.6 (0.183) (0.344) = 3.63 for 0^ 

= 57.6 (0.183)(0.814) = 8.58 for 0g 


(28) 


Group 3: 


= 6.50 - (1.45 + 3.63) = 1.42 


(29) 


AW3 = 


0.329 Yg 


1.42 


57.6(0.187) 


= 0.132 




Beneath diagonal line 
in 0™ column 

m3 


57.6 (0.132) (0.187) = 1.42 for 64 

57.6 (0.132) (0.351) = 2.66 for 65 


57.6 (0.132) (0.704) = 5.35 for 0g 
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Group 4: 



= 9.78 - (1.79 + 5.32 + 2.66) = 0.01 


''0 


= 0.00095 


and so forth. This process is continued until a value of Av is deter- 
mined for each group. 

The resulting "block diagram is shown "by the straight solid lines in 
figure 26. A solid curve is faired through this block diagram. For a 
known value of . (p, figure 15 may be used to determine K and hence a., 
droplet radius, as a function of 0jjj. The droplet size scale thus deter- 
mined is shown along with the 0jjj scale in figure 26. The values of w 

as a function of droplet radius read from the solid curve of figure 26 
are identical to those for one of the three 2-inch-cylinder curves in 
figure 21(c). 

The significance of the increments selected for row 1 in the 

preceding table is illustrated by the dotted curve in figure 26, which 
results from dividing the Wp curve at 0 = 83.7, 75, 65, 55, 45, 35, 

25, 15, 5, and 0 degrees. The dashed and the solid curves are almost 
identical. 
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Figure 2. - Hypothetical dimensionless distribution of droplet sizes. 
Langmuir D-type distribution from reference 2. 
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Cumulative fraction of total liquid-water content, n 
(Decreasing droplet size -► ) 

(a) Block- type droplet distribution. 

Figure 3. - Typical representation of weighted total collection 
efficiency Ej^ as summation of contributions by each droplet 
size in distribution. Cross-hatched area equals 
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Cumulative fraction of total liquid-water content, n 
(Decreasing droplet size -*■ ) 

(b) Faired-type droplet distribution. 


Figure 3. - Concluded. Typical representation of weighted total 
collection efficiency Em as summation of contributions by _ 
each droplet size in distribution. Cross-hatched area equals E 
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0 .2 -4 *6 -o 1 

Cumulative fraction of total liquid water content, n 
(Decreasing droplet size -*• ) 

Figure 4. - Typical representation of weighted local collection 
efficiency P as sunnnation of contributions by each droplet- 
size group in distribution for one particular surface location 
Sq on two-dimensional body for_b lock- type droplet distribution 
Cross-hatched area represents P at Sq- 
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Cumulative fraction of total liquid -water content, n 
(Decreasing droplet size -► ) 

Figure 5. Typical variation of weighted local collection effi 
clency p (represented by area under a particular s^ curve 
as function of surface location on two-dimensional body for 
block-type droplet distribution. 
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Increasing surface distance from reference point 

Figure 6. - Typical variation of weighted local collection eff 
^ as function of surface location on two-dimensional body f 
block-type droplet distribution. 
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Cumulative fraction of total liquid-water content, n 
(Decreasing droplet size -►) 

Figure 7. Typical variation of weighted local collection effi- 
ciency p (represented hy area under particular 8^ curve) as 
function of surface location on two-dimensional body for faired- 
type droplet distribution. 
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Increasing surface distance from reference point 


Figure 8. - Typical variation of weighted local collection 
efficiency P" as function of surface location on two- 
dimensional body for faired-type droplet distribution. 
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Figure 9. - Typical installation of cylinder (4- in. diam.), 
with blotter attached, in 6- by 9-foot test section of 
icing tunnel (looking downstream). 
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Figure 10. - Installation of grid used for mounting ribbons to measure 
uniformity of spray cloud in 6- by 9- foot test section of icing 
tunnel (looking downstream). 
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Figure 11. - Cross section of air- water atomizing-nozzle assembly. 
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Figure 13. - Trajectory starting ordinates as function of angle of impingement 
for cylinder (ref. l) . 
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Figure 14. - Local collection efficiency for cylinder 
as function of central angle (ref. 1) . 
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Figure 15. - Limit of impingement on cylinder (ref. l) . 
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Figure 18. - Block- and falred-type droplet size dlS' 
tributlons as specified by limit of impingement 0, 
on cylinder and cumulative liquid-water content w 
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(a) Air-water pressure ra 
120 psl water) . 


Figure 20. - Experimental 
three cylinder sizes, j 
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0 20 40 60 80 100 

Central angle of cylinder, 0, deg 


(b) Air -water pressure ratio, 0.6 (60 psi air to 
100 psi. water) . 

Figure 20. - Continued. Experimental local water impingement 
rates for three cylinder sizes. Airspeed, 175 miles per 
hour. 
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0 20 40 60 80 100 

Central angle of cylinder, 0, deg 


(c) Air-water preseure ratio, 0.8 (80 pel air to 100 psi water) . 

Figure 20. - Concluded. Ejqperimental local water In^iingement 
rates for three cylinder sizes. Airspeed, 175 miles per hour. 
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Figure 21. - Continued. Droplet size distribution obtained from experimen- 
tal impingement rates on cylinders (based on theoretical trajectories 






(c) Air -water pressure ratio, 0,8 (80 psi air and 100 psl water). ‘ 

Figure 21. - Concluded. Droplet size distribution obtained from experimental Impinge 
ment rates on cylinders (based on theoretical trajectories from ref. l) . 
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Figure 22. - Droplet size distribution obtained from experimental ln5)lngement rates on cylinders 
(based on theoretical trajectories from ref. 2) . 
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Droplet radius^ a, microns 

(c) Air-water pressure ratio, 0.8 (80 psl air to 100 psl water). 

Figure 22. - Concluded. Droplet size distribution obtained from experimental impingement 
rates on cylinders (based on theoretical trajectories from ref. 2) . 
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Droplet radius, a, microns 


Figure 23. - Final evaluation of droplet size distribution for three spray 
conditions based on experimental Impingement rates of cylinders and an 
aspirating device (based on theoretical trajectories from ref. l) . 
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Figure 26. - Block- and falred-type droplet size distribution as specified by limit of Impingement (also 

droplet radius a) and cumulative liquid-water content w from numerical exanple in appendix C.- 




